Evidence for AGN Feedback in Galaxy Clusters and Groups 
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, The current generation of flagship X-ray missions, Chandra and XMM-Newton, has changed our 

■ understanding of the so-called "cool core" galaxy clusters and groups. Instead of the initial idea 
Qj^' that the thermal gas is cooling and flowing toward the center, the new picture envisages a complex 

dynamical evolution of the intra-cluster medium (ICM) regulated by the radiative cooling and the 
nongravitational heating from the active galactic nucleus (AGN). Understanding the physics of the 
'/^ ' hot gas and its interplay with the relativistic plasma ejected by the AGN is key for understanding 

the growth and evolution of galaxies and their central black holes, the history of star formation, 
and the formation of large-scale structures. It has thus become clear that the feedback from the 
central black hole must be taken into account in any model of galaxy evolution. In this paper, we 
draw a qualitative picture of the current knowledge of the effects of the AGN feedback on the ICM 
, by summarizing the recent results in this field. 
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' The physics of the intra-cluster medium (ICM) of clusters and groups of galaxies is complex. The current 

I generation of X-ray satellites, Chandra and XMM-Newton, has shown indeed that it is regulated by yet poorly 

understood non-gravitational processes beyond simple gravity, gasdynamics and radiative cooling usually 
considered in the standard cold dark matter cosmological scenario [225l |. In particular, an important discovery 
from high-resolution X-ray observations was that the amount of thermal gas radiatively cooli ng to low 
' temperatures is much less than what is predicted by the standard "cooling flow" model (see [s^ [Sll, Il72l Il73l | 

ff^ , and references therein), thus radically changing our understanding of the so-called "cool core" systems. The 

CO ' implication is that the central gas must experience some kind of heating due to a feedback mechanism that 

CO I prevents cool cores from establishing cooling flows at the rates predicted by earlier, low-resolution X-ray 

^\ • observations. Establishing the source of this heating, and understanding when and how it takes place, has 

' become a major topic of study in extragalactic astrophysics. 

Based on observational evidence and theoretical modelling, the primary source of feedback has been 
identifled in the outbursts and accompanying ene rgy injection, likely intermittent, from the active galactic 
nucleus (AGN) of the dominant cD galaxies (e.g., |145j and references therein), which host the most massive 
black holes in the local Universe. AGNs manifest as central radio sources, which are commonly observed 
. , in cool core clusters Most of the cool core systems have highly disturbed X-ray morphologies, and 

■ radio observations clearly show that AGN jets are the cause of many of the structures revealed by the X- 
ray telescopes. Such surface brightness features, including apparent depressions or "cavities" in the X-ray 
images and sharp density discontinuities interpreted as shocks, indicate a strong interaction between the 
central AGN and the intra-cluster medium (ICM). The incidence and variety of cavities, shocks, and ripples 
observed both in the radio and in X-r ays i n the hot ICM provides direct evidence of the widespread presence 
of AGN-driven phenomena (e.g., (TtI , llOOf for sample studies of clusters and groups, respectively). 

Such AGN feedback has a wide range of impacts, from the formation of galaxies, through to the explanation 
of the observed relation between the black hole mass and the bulge velocity dispersion (which indicates a 
caus al co nnection or feedback mechanism between the formation of bulges and their central black holes, 
e.g., [134 ), to the regulation of cool cores which explains why cooling and star formation still proceeds at a 
reduced rate. On the other hand, the details of how the feedback loop operates are still unknown. Feedback 
is also required to suppress the overproduction of massi ve g alaxies predicted by dark-matter-only simulations 
and to break the self-similarity of cluste rs [ e.g.. [lOl IT^, l57j[. The nature of this feedback is therefore vital to 
our understanding of galaxy evolution ('49] and references therein). 

After a brief discussion of the importance of galaxy clusters and their scaling relations for the study of 
cosmic evolution (|jlT|, we overview the role of AGN feedback in structure formation ( mil\i and the basic 
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properties of clusters in X-rays eraphasizing the hot intra-cluster medium fi jlVp . We then focus on the 
observational evidence of AGN feedback in action in galaxy clusters and groups (SjVj and finally give our 
conclusions ( Wip . The present paper does not intend to be a comprehensive review, but aims instead at 
drawing a qualitative picture of the impact of the A GN feedback on the ICM addressed particula rly t o the 



novices in this field. General revie ws of clusters from an X-ray perspective were given by Sarazin |195| and 
more recently by Mushotzky (l58j | and Arnaud [3]. A review of clusters as cosmological probes was given 
by Voit [219], and cold fronts and shocks associated with cluster mergers were reviewed by Markevitch and 
Vikhlinin [133] ■ An exhaustive r evie w of the issues of AGN heating in the hot atmospheres was recently 
given by McNamara and Nulsen jl45j ]. 

Throughout the paper we assume a cosmology with Hq — 70 km s^^ Mpc^^ and il^/ = 1 — f^A = 0.3, 
where not specified otherwise. 



II. CLUSTERS OF GALAXIES AS COSMOLOGICAL PROBE 



The existence of clusters of galaxies and of other cosmic structures demonstrates that the Universe is not 
perfectly homogeneous. The matter density of the primordial Universe must have been slightly inhomoge- 
neous, with overdense perturbations which deviate from the mean density. In the so-called "Concordance 
Model" largely accepted today as the standard cold dark matter cosmological scenario, cosmic structures like 
galaxies and clusters of galaxies originate from the gravitational instability of these primordial density fluc- 
tuations. The formation of structures from perturbations in the density distribution of cold dark matter is a 
hierarchical process. Small subclumps of matter are the first to deviate from the Hubble flow, to collapse and 
to experience gravitational relaxation because the density perturbations have larger amplitudes on smaller 
mass scales. These sma ll ob jects then undergo a merging process to form larger and larger structures, up to 
the clusters of galaxies (225| . 

Galaxy clusters trace the high-density tail of the primordial field of dark matter density perturbations, 
and their num erica l density as a function of redshift, z, is highly sensitive to the specific cosmological 
scenario (e.g., |l89 l| and references therein). Therefore, if one builds the so-called "cluster mass function" 
nMiM), i.e. the number density of clusters with mass greater than M in a comoving volume element, and 
determines its evolution with redshift, it is possible to constrain the main cosmological parameters from 
the comparison between the observations and the theory predictions. A complete understanding of the 
details of the process of hierarchical merging would require accurate numerical simulations. However, many 
fundamental aspects can be illustrated by spherically-symmetric, simplified models of cluster formation. 
In particular, the combination of spherical top-hat collapse models with the growth function for linear 
perturbations (assumed to be gaussian) has led to a variety of semi-analytical methods t o exp ress the cluster 
mass function in terms of cosmological parameters (the seminal work in this field is by [179j). 

Therefore, the comparison between the theoretical mass function and the mass function determined from 
observations allows one to constrain the main cosmological parameters, although with a degeneracy between 
the matter density parameter, VIm, and the power spectrum normalization of the perturbations within a 
comoving sphere of radius 8 Mpc, cts- Such a degeneracy can be broken by studying the redshift 
evolution of the mass function, which is highly sensitive to 57^/, by taking into account the evolution of 
the observables (see the reviews by [2^, Il89l . |219| and references therein) . The accuracy of the cosmological 
parameter measurements is currently limited by uncertainties in the relations between cluster masses and 
the observable properties that trace these masses, such as luminosity or temperature. In order to measure 
the mass function from a large sample of clusters is indeed necessary to link the mass to these quantities 
which are easily observable. 

In this context are very useful the so-called "self-similar scaling relations", derived naturally by considering 
that the cosmological structures originate from scale-free density perturbations and that the thermodynam- 
ical properties of the ICM are determined by scale- free gravity only jl22l ]. Under these assumptions, galaxy 
clusters of different masses may be considered as a scaled version of each other. The density of each dark 
matter halo, pdm, is proportional to the critical density of the Universe at the cluster's redshift, Pc,z, through 
the so-called "overdensity" A = pdm/Pc,z, where pc,z = /SttG, and the expression for the Hubble con- 
stant at redshift 2 in a flat ACDM Universe is [e.g.,[38S]: = Hox/fimC^ + z)^ + fl^ ^ Ho E{z). Therefore, 
all clusters should have the same properties when rescaled by A. 

If we define the mass M as the mass Ma inside the radius i?A at a given overdensity A, we obtain: 
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Ma oc Pc^z ■ A • cx pc.o ■ E{zY ■ A ■ i?^, and thus we get to the M ~ R relation in the form 

i?cx Ml/3. -2/3 (1) 

During cluster formation, the gravitational collapse of the diffuse gas in the potential well of the dark matter 
halo heats the gas itself at the virial temperature of the potential well that confines it: T^--^ ~ GM fim-p/kR ^ 
lO^K where M is the total mass, k is the Boltzmann constant, fi ^ 0.6 is the mean molecular weight and 
R is the virial radius. The gas is thus heated to X-ray emitting temperatures and becomes a plasma in 
hydrostatic equilibrium whose emissivity is proportional to the square of its density (see i jlVp . The virial 
temperature of an isothermal sphere of mass AI is: kT cx M/R oc M^/^ • E{z)'^^'^, leading to the M — T 
relation in the form 

M (X r3/2 . E{z)-^ (2) 

From these relations it is possible to derive the relation between temperature and luminosity emitted by 
the hot gas through bremsstrahlung emission: Lx oc ■ A-V, where p is the average gas density and A is the 
cooling function, that in the bremsstrahlung regime is oc T^/^ (see ijlVp . Assuming that the gas distribution 
traces the dark matter distribution, p oc poM oa Pc,z, we obtain: Lx oc p ■ T^/"^ ■ M (x po ■ E{z)'^ ■ T^/^ ■ M oc 
£:(z)2 . Ti/2 • T3/2 . E{z)-^ thus deriving the L - T relation in the form: 

Lx^T^- E{z) (3) 

By combining the M — T relation (Eq. [J) with the L — T relation (Eq. [3]), we can finally derive the M — L 
relation that links the mass directly to the observable luminosity: Lx oc [A-'l'^/'^E{z)'^^'^]'^E{z), finding 

MocL^/^-S(z)-^/4 (4) 

In principle, once calibrated with simulations and/or observations, these scaling relations provide a method 
to link the mass of clusters to observables under the assumption that the process of structure formation is 
guided by gravity alone. On the other hand, deviations from these relations testify the presence of physical 
processes more complex than gravitational dynamics only, which modify the thermodynamical properties 
of the diffuse baryons and therefore the relations between observables and cluster masses. In particular, a 
number of observational measurements seems to indicate that the L — T relation is s teeper that predicted 
by self-similar models, and is in the form L oc T^-^-^ (e.g., 0, [t^, [ill [iM Ell [izS . see also WF\i . This 
observed breaking of the scaling relation has been ascribed to the presence of some excess entropy in the 
gas due to primordial nongravitational heating before the cluster virialization [79, 123], and is one of the 
strongest evidence for nongravitational processes acting in the ICM. 

The main source of uncertainty in the determination of cosmological parameters from studies of cluster 
samples arises then from the uncertainty in the normalization, shape and evolution of the relationships that 
relate the cluster masses to the observables. In order to understand better such relations it is essential to 
investigate how the structure formation and AGN feedback affect the evolution of what we can observe, i.e. 
the baryons in clusters. 



III. ROLE OF AGN FEEDBACK IN GALAXY EVOLUTION 



O ne o f the main problems of the current cosmological model is why so few baryons have formed stars 
[5^ |224| |. Numerical simulations of cosmological structure formation that include the hydrodynamics of 
baryons and the radiative cooling processes predict that ^20% of the baryons should have condensed into 
galaxies, but only ^10% have been observed in the form of stars [e.g.,Q. In particular, simulations that 
include only gravitational heating predict an excessive cooling of baryons that results in a population of 
galaxies which are too massive and too bright with respect to the ones observed, thus failing to reproduce 
the truncation of the high-luminosity end of the galaxy luminosity function JA,, .199). 

Instead of residing in the cD galaxies as predicted by simulations, most baryons are observed in the 
hot ICM. This problem may find a solution in the nongravitational heating supplied by supernovae (SN) 
and active galactic nuclei (AGN). Supernovae are essential in the process of enrichment of the ICM to the 
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metallicity level observed [29l llSOj , and from the heavy-element abundances in clusters it is estimated that 
during a cluster history they supply a total amount of energy of the order of 0.3-1 keV/particle [e.g..l90l.[l76t. 
This is not enough energy to quench cooling in massive galaxies , as the energy input required to explain 
the mass-observable relations is ~ 1 — 2 keV/particle |.211 , 221 , .229.] . Energetically, AGN heating appears to 
be the most likely mechanism to severely reduce the supp l y of gas from the hot ICM in massive galaxies and 
to explain the observed entropy profiles [ll [gi, [HI [Hi H^l • AGNs arc powered by accretion of material 
onto a black hole (BH) , which is located at the center of each stellar spheroid (both bulges within spirals and 
ellipticals). Matter falling onto a black hole releases an energy of the order of Ebh = eMc^, where e « 0.1 
is the efficiency. For supermassive black holes (SMBH) of masses ^ IO^Mq, the amount of energy released 
during their formation and growth is of the order of Ebh ~ 2 x 10^^ erg s^^. Even a tiny fraction (^1%) of 
the energy released within the bulge could heat and blow away its entire gas content in small systems and 
prevent cooling, thus explaining the lack of star formation in bulges. An extraordinary discovery obtained 
recently in astrophysics is the correlation between the mass of the central black hole (Mrh) and the velocity 
dispersion (cr) of the galaxy's bulge, used to estimate the mass of the bulge itself [134l|. This "Magorrian 
relation" A^bh — cr suggests that the large-scale properties of the galaxy and the small-scale properties of 
the black hole are related. In particular, each massi ve g alaxy seems to host a central black hole, whose 
mass is ~ 0.1 — 0.2% of the bulge stellar mass [sol. [osl [2T2I . Such a correlation may arise from the fact that 
the central black hole is able to regulate the amount of gas available for star formation in the galaxy. The 
formation of black holes and the formation of bulges are closely linked. Therefore supermassive black holes 
can have a profound influence on the formation and evolution of galaxies. The physical process regulating 
these phenomena has been called "feedback", and the understanding of how it acts in detail is one of the 
main open issues in extragalactic astrophysics. 

Clusters of galaxies are the only locations in the Universe where we can find an almost complete census of 
the intergalactic baryons and a very good description of their thermodynamical status and of their enrichment 
in heavy elements. Therefore, X-ray observations of the ICM can provide us with new important insights 
into the processes of cooling and feedback which regulate galaxy formation. 



IV. CLUSTERS OF GALAXIES IN X-RAYS AND THERMAL ICM 



Clusters of galaxies are the largest virialized structures in the Universe, with typical sizes of 2 — 4 
Mpc and total gravitational masses of ~ 10^"* — IO^^Mq. They are luminous X-ray sources, with typical 
luminosities ranging from a fewx 10^'^ — 10^^ erg s~^. As first suggested by Felten et al. ^1], the X-rays from 
clusters are primarily thermal bremsstrahlung emission from the diffuse ICM which fills the deep potential 
wells and is heated to temperatures of ~ 10^ K (where kT — 1 keV for T = 1.16 x 10^ K) during the process 
of cluster formation. 



A. Physical Properties of Hot Diffuse Plasma 



The simple assumptions which are generally made in the study of the ionizat ion s tate and X-ray line and 
continuum emission from a low density, hot plasma are briefly reviewed below '195*]. 

1) The time scale for elastic Coulomb collisions between particles in the plasma is much shorter than the 
age or cooling time of the plasma, therefore the free particles are assumed to have a Maxwell-Boltzmann 
distribution at the temperature T. This follows from considerations on the mean free paths of particles in a 
plasma without a magnetic held. The mean free pa th A c for an electron to suffer an energy exchange with 
another electron via Coulomb collisions is given by [205j : 

_ s^Hkr^r ... 



where To is the electron temperature, ric is the electron number density and A is the ratio of largest to smallest 
impact parameters for the collisions (In A w 38). Eq. [5]assumes that the electrons have a Maxwellian velocity 
distribution at the electron temperature Te. However, it can be demonstrated that if a homogeneous plasma 
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is created in a state in which the particle distribution is non-Maxwelhan, elastic colli sions will cause it to 
relax to a Maxwellian distribution on a time scale determined by the mean free paths (205l l206j | . Electrons 
will achieve this equilibrium on a time scale given roughly by toq(e, e) = Ac/('yc)rmsi where ('yo)ims is the rms 
electron velocity = (S/cTo/rrio)^^^: 

The time scale for Coulomb collisions between protons to bring them into kinetic equilibrium is about 
tcq{p,p) ~ (TOp/mc)^/^icq(e, e), roughly 43 times longer than that for electrons. After this time, the electrons 
and ions (generally assumed to be protons) would each have Maxwellian distribution, but generally at 
different temperatures, respectively To and Tj. The time scale for the electrons and ions to reach equipartition 
To = Ti is <oq(p, e) ~ (mp/mo)ioq(e, e), and for typical values of the ICM temperature and density is teq{p, e) 
<6 X 10* yr. Since this is shorter than the age of the clusters or their large-scale cooling time (although it 
is comparable to or longer than the cooling time in the cores of clusters, see Eq. \T\\ below) . the intra-clustcr 
plasma can generally be characterized by a single kinetic temperature T = = Tj, which determines the 
rates of all excitation and ionization processes. It is important to note that the mean free paths, i.e.: 

Ae = Ai«23kpcf^vVf (7) 

are generally much shorter than the length scales of interest in clusters (w 1 Mpc), and therefore the ICM 
can be treated as a coUisional fluid, satisfying the hydrodynamic equations. 

2) At these low densities, collisional excitation and de-excitation processes are much slower than radiative 
decays, therefore any ionization or excitation process is assumed to be initiated from the ground state of an 
ion. Three (or more) body collisional processes are ignored because of the low density. 

3) Stimulated radiative transitions are not important, since the radiation field in the ICM is sufficiently 
dilute. 

4) At these low densities, the gas is optically thin and the transport of the radiation field can therefore be 
ignored. 

Under these conditions, ionization and emission result primarily from collisions of ions with electrons, and 
collisions with other ions can be ignored. The time scales for ionization and recombination are generally 
considerably shorter than the age of the cluster or any relevant hydrodynamic time scale, therefore the plasma 
is assumed to be in ionization equilibrium [e.g., |203|. The equilibrium ionization state of a diffuse plasma 
depends only on the electron temperature: since in nearly all astrophysical plasmas most of the electrons 
originate in hydrogen and helium atoms, and these are fully ionized under the conditions considered here, 
the ICM is generally treated as a fully ionized plasma. 

By indicating with X, y, Z the mass fraction of hydrogen, helium, and heavier elements, respectively, 
the corresponding atom number densities can be written in the form: njj = rip = pX/rrip, njie = pi^/4mp = 
UpY/AX, riz — pZ/Anip — UpZ/AX, where p is the gas density, nip the proton mass and A the mean atomic 
mass number (i.e., the number of nucleons) of heavier elements. Assuming that the gas pressure p = nkT 
is contributed only by electrons and protons, thus neglecting nuclei (n = rio + Up), it is possible to derive 
the electron density rie in terms of the proton density rip. From the expression of the number of particles 
contributing to the pressure, n — 2nH + Snfjo + ^An^, one obtains: 

/ lY 1Z\ 

which for solar abundances {X — 0.71, Y ~ 0.265, Z ~ 0.025) leads to Uc ^ 1.2np. It is also possible to 
calculate the mean molecular weight in amu, p,, such that the total number density of particles (electrons, 
protons and ions) is n = p/ pnip. From the expression n = 2nH + 3nHo + (^A + l)nz, in the approximation 
A:S> 1, one obtains: 
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M-(2x + ^r + iz^ (9) 

which for solar abundances leads to /i ^ 0.6. 

B. X~ray Emission from the ICM 

The X-ray continuum emission from a hot diffuse plasma, such as the ICM, is due primarily to two 
processes: thermal bremsstrahlung (free-free emission) and recombination (free-bound) emission. Processes 
that contribute to X-ray line emission (bound-bound radiation) from a diffuse plasma include coUisional 
excitation of valence or inner shell electrons, radiative and dielectric recombination, inner shell coUisional 
ionization and radiative cascades following any of these processes. 

At the high temperatures typical of clusters (in particular at kT ^2.5 keV), thermal bremsstrahlung is 
the predominant X-ray emission process. The bremsstrahlung emissivity at a frequency v (defined as the 
emitted energy per unit time, freq uency and volume) of a plasma with temperature T, electron density 
and ion density rii is given by [e.g.. Il94l |: 

Jbr(i^, T) = 6.8 X 10-38z2nenir-i/2e-'*''/'=^5(z/, T) (10) 

where the Gaunt factor g{T)^ which corrects for quantum mechanical effects and for the effect of distant 
collisions, is a slowly varying function of the parameters 124, 125]. If the ICM is mainly at a single 
temperature, then Eq. [10] indicates that the X-ray spectrum should be close to an exponential of the 
frequency, as is generally observed. 

The total power per unit volume emitted by thermal bremsstrahlung is obtained by integrating Eq. 1101 
over frequency, obtaining: 

Jb^(T) = 1.4 X lO^^"^ ncniT^/^Z^g{T) (11) 

where g(r) is a frequency average o{g{i>,T) and is in the range 1.1 to 1.5 (choosing a value of 1.2 will give 
an accuracy in the estimate of Jbr{T) to within about 20%, 1941). For solar abundances, the emission is 
primarily from hydrogen and helium. 

C ompi lations of the different emissivities for X-ray lines and continua can be found in the literature [e.g., 

Il5ll Il85j | . Early detailed calculatio ns o f the X-ray spectra predicted by different models of the ICM have 
been given by Sarazin and Bahcall Jl96!] and Bahcall and Sarazin [3|. In these models most of the X-ray 
emission is thermal bremsstrahlung continuum, and the strongest lines (highest equivalent width) are in the 
7 keV iron line complex. 

The frequency-integrated total emissivity at a temperature T can be written as: 

Jx(r) = A(T)nenp erg s^^ cm-^ (12) 

where A(T) is the cooling function, which depends upon the mechanism of the emission and can be represented 
as: 

A(r) = /T" (13) 

where —0.6 < a <0.55; for thermal bremsstrahlung it is Z ~ 2.5 x 10"^"^ and a = 1/2 Il44ll. The general 
behaviour of the cooling function was calculated and discussed by Sutherland and Dopita |209l | . 

The projection on the sky of the plasma emissivity gives the X-ray surface brightness: in order to constrain 
the physical parameters of the ICM, the observed surface brightness can be either geometrically deprojected 
or, more simply, fitted with a model obtained from an assumed distribution of the gas density. 
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C. Hydrostatic Models for ICM Distribution 

From the expression for the sound speed Cg^ — ^kT / ^rrip, where 7 — 5/3 for a monatomic gas, the time 
required for a sound wave in the ICM to cross a cluster is given by: 

_ 1 

-,8 / \ I ^ 



ts ~ 6.6 X 10"* — yr (14) 

where D is the cluster diameter. Since this time is short compared to the likely age of a cluster (in first 
approximation assumed to be ^ 10^° yr), the gas is generally thought to be in hydrostatic equilibrium in 
the gravitational potential of the cluster: Vp — —pVcj), where p — pkT / ^m-p is the gas pressure, p is the gas 
density and is the gravitational potential of the cluster. Under the assumptions that the ICM is locally 
homogeneous and the cluster is spherically symmetric, the hydrostatic equilibrium equation reduces to: 

Idp ^ ^ GMjr) 
p dr dr 

where r is the radius from the cluster center and M{r) is the total cluster mass within r. If the gas self 
gravity is ignored, then the distribution of the ICM is determined by the cluster potential (j){r) and the 
temperature distribution of the gas T{r), and Eq. [T5]is a linear equation for the logarithm of the gas density. 
Under these assumptions, the gravitational mass Aftot of a galaxy cluster can be written as: 



, , / X kT r 
Mtot(<r)--- 

G/iTOp 



din n d In T 

H 

d In r d In r 



(16) 



This expression is commonly used to estimate the gravitational mass of galaxy clusters and groups from 
X-r ay observ ations, through the measurements of the radial profiles of temperature and density [e.g., ItsI 
[96l . llOSl |217| . However, we note that Eq. [TBI neglects the contribution of possible additional, non-thermal 
pressure that, if present, should be included in the estimate of the total mass. In particular, recent results 
from numerical simulations indicate that the total mass of simulated clusters estimated through the X-ray 
approach is lower that the true one due to gas b ulk motio ns (i.e. deviation from the hydrostatic equilibrium) 
and the complex thermal structure of the gas [l59l Il84j . P ossib le observational biases in the derivation of 
X-ray masses are also discussed in Piffaretti and Valdarnini jl75j |. 



The 13 -Model 



Cavaliere and Fusco-Femiano (48| studied the X-ray emission by the hot plasma in galaxy clusters and 
developed a hydrostatic model based on the assumption that the gas and the galaxies are in equilibrium 
in the same gravitational potential <j) (see e.g., ^95,j, and ^ for a recent commentary on this model). By 
further assu ming that the galaxy distribution is well described by King's approximation to the isothermal 
sphere jl27l |. the expression for the ICM distribution may be written as [48j : 



Pir) = Po 



1 



(17) 



and the surface brightness profile observed at a projected radius 6, I{b), is in the form |48| : 



I{h) = lo 



1 



The parameter /3 is defined as 



13^ 



kT/p,mp 



(18) 



(19) 
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where is the Une-of-sight velocity dispersion, and represents the ratio of specific kinetic energies of galaxies 
and gas. 

This self-consistent isothermal model, called the "(3 -model", is widely used in the X-ray astronomy to 
parametrize the gas density profile in clusters of galaxies by fitting their surface brightness profile. One of 
the advantages of using a /3-model to parameterize the observed X-ray surface brightness is that the total 
mass profiles can be recovered analytically and expressed by a simple formula: 



Mtot{< r) 



GflTTlp 



3(3rT dT 

J.2 _|_ J.2 ^j. 



(20) 



Eq. [T7] states that the gas density rises towards the cluster center. Since the bremsstrahlung and line 
emission depend on the square of the gas density (Eq. 1121) , in the central regions of clusters the loss of energy 
by X-ray emission represents an important process for the thermal particles in the ICM. In particular, if the 
gas density reaches high enough values, large amounts of gas cool and flow into the centers of clusters, forming 
the so-called cooling flows. In cooling flow clusters, the single /3-model is found to be a poor description of 
the entire surface brightness profile: a fit to the outer regions shows a strong excess in the center as compared 
to the model (see ^IVPp . Conversely, the centrally peaked emission is a strong indication of a cooling flow 
in relaxed cluster. 



D. Cool Cores and Cooling Flow Problem 

The X-rays emitted from clusters of galaxies represent a loss of energy of the ICM. The resultant cooling 
time is calculated as the time taken for the gas to radiate its enthalpy per unit volume using the 
instantaneous cooling rate at any temperature: 

Hy _ J kT 

where: 7 = 5/3 is the adiabatic index; /i w 0.6 is the molecular weight; X « 0.71 is the hydrogen mass 
fraction; and A(r) is the cooling function. In the central region, the cooling rate of the ICM is sufficiently 
high that the particles lose their energy via radiation, as inferred from X-ray images of the cores of many 
clusters which show strongly peaked surface brightness distributions. The density of the gas then rises to 
maintain the pressure required to support the weight of the overlying gas in the rest of the cluster, causing a 
slow subsonic inflow of material towards the cluster center. This qualitative picture describes the physics of 
the process known as a cooling flow (see [si'] for a review of the standard model, and ^76^ for a quantitative 
description of the evolution of cooling flows). The steady cooling flow is confined within the region in which 
tcooi is less than the time for which the system has been relaxed (thus allowing the establishment of a cooling 
flow). This cooling region is delimited by the so-called cooling radius Tcoob which is usually defined as the 
radius at which <cooi is equal to the look-back time to z = 1, i.e. ^ 7.7 x 10^ yr in the concordance cosmology. 
The fraction of clusters with a central surface excess with respect to a /3-model, the so-called cool cores, is 
large. Cool core clusters are about 90% of X-ray-selected clusters with total mass Mtot <1O^^M0, and about 
50% of X-ray-selected clusters with total mass Mtot ^lO-^^'M© [49|. Cool cores are also characterized by 
strong enhancements in th e central abundance [e.g.. IgiI. l62j and declining temperature profiles toward the 
central region [e.g.. [Tl. [216j . 

In the standard model, the "magnitude" of a cooling flow is measured from the amount of matter which 
crosses rcooi and flows towards the center, that is M, the mass inflow rate. The mass inflow rate, due to 
cooling, can be estimated from the X-ray images by using the luminosity Lcooi associated with the cooling 
region and assuming that it is all due to the radiation of the total thermal energy of the gas plus the pdV 
work done on the gas as it enters re', icooi = dE/dt, where dE = dEth + pdV = (7/7 — l)pdV, and 
pdV = {p kT dV) / (finip) — {dM kT) / {finip) , with 7 = 5/3. By substituting one obtains the expression for 

-^cool- 



(22) 
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where T is the temperature of the gas at TcooI- Lnoo\ ranges from ~ 10"*^ to > 10'*'* erg s~l and generally 
represents ~ 10% of the total cluster lummosity [Slj. Value of M ~ IOOMq yr~* are fairly typical for cluster 
cooling flows. 

However, the current generation of X-ray satellites, Chandra and XMM-Newton, has radically changed our 
understanding of cooling flow systems. Albeit confirming the existence of short cooling times, high densities 
and low temperatures in the cluster cores, the arrival of high-resolution X-ray spectral data has shown the 
absence or weakness of the soft X-ray line Fe XVII, indicating that the amount of gas coo l ing radiatively 
below about one third of its original temperature is ten times less than expected [e.g-. fssllllSl . 11731 . The lack 
of evidence for central gas cooling to very low temperatures at the rates predicted in the hot atmospheres 
of galaxy clusters and groups re presents an open question which is often referred to as the so called 'cooling 
flow problem' (see psi . Il45l Il72j | for reviews). 

Historically, two main approaches were adopted to solve this problem. As the gas radiates but does not 
appear to cool, either the normal signatures of radiative cooling below 1 — 2 keV are somehow suppressed, or 
there must be an energy-injection mechanism into the ICM w hich compensates cooling. Differei it possib ilities 
considered in the former hypothesis include absorption [13, Il73j ) , inhomogeneous metallicity [SJ, Il55( , and 
the emerging of the missing X-ray lum inos ity in other bands, like ultraviolet, optical and infrared due 
to mixing with cooler gas/dust [SJ, [H, Il39l |. Proposed heating mechanisms in the c ontext of the latte r 
approach include e.g., processes associated with relativistic AGN ou tflows Isil . 33, [HI, Il2ll . Il90l Il9ll l210j |. 
electron thermal conduction from the outer regions of clusters [86l.l213L[2lf l23a|. continuous subcluster 
merging (l38[. contribution of the gravitational potential of the cluster core |82| , feedback from intra-cluster 
supernovae [65|. Among all these, feedback by the central AGN appears to be the most promising solution. 



V. X-RAY CAVITIES AND SHOCKS: AGN FEEDBACK IN ACTION 



It was already known in the early 90s that central dominant (cD) galaxies of cool core clusters have a high 
incidence of radio activity, showing the presence of central FR-I radiogalaxies in 70% of the cases 
Their behaviour differs from that of quasar: in many low-accretion-rate AGNs almost all the released energy 
is channelled into jets because the density of the gas surrounding the black hole is not high enough for an 
efficient radiation [e.g., H^]. In fact, the importance of these objects has been underestimated for a long 
time due to their poor optical luminosity. The importance of radio galaxies in cool cores began to emerge 
after the discovery, with the X-ray satellite ROSAT, of deficits in the X-ray emission of the Perseus and 
Cygnus A clusters which are spatially coincident with regions of enhanced synchrotron emission [U, l42j . 
With the advent of the new high-resolution X-ray observations performed with the current generations 
of X-ray telescopes, Chandra and XMM-Newton, it became clear that the central radio sources have a 
profound, persistent effect on the ICM. In particular, Chandra images, which are obtained at the superb 
angular resolution of '--^l", showed that the Perseus and Cygnus A clusters are not isolated cases - indeed 
the central hot gas in many cool core systems is not smoothly distributed, but is instead cavitated on scales 
often approximately coincident with lobes of extended radio emission. These observations also reveal highly 
disturbed structures in the cores of many clusters, including shocks, ripples and sharp density discontinuities. 
The comparison with radio images obtained at similar angular resolution has revealed that AGN jets are the 
cause of many of these disturbances. The most typical configuration is for jets from the central dominant 
elliptical of a cluster to extend outwards in a bipolar flow, inflating lobes of radio-emitting plasma. These 
lobes push aside the X-ray emitting gas of the cluster atmosphere, thus excavating depressions in the ICM 
which are detectable as apparent 'cavities' in the X-ray images. The brightness depressions observed in 
X-rays, which are mostly found in pairs, are ^ 20 — 40% below the level of the surroun ding gas, consistently 
with the expected decrement along an empty bubble immersed in a /3- model atmosphere [I9l[20ll83l . ll46lll62l |. 
The cavities are often surrounde d by bright shells, or rims, which are typically found to be cooler than the 
ambient medium jl^, [6^ [H, Il46{ . This is likely due to the compression of the central, low-entropy gas 
into the bright shell during the cavity rising and expansion into the hot atmosphere. 

X-ray cavities are present in >70% of cool core clusters [tO], but this fraction could be underestimated 
due to the limitation of cavity detectability [l^ ■ Identifying radio galaxies as a primary source of feedback 
in the hot atmospheres of galaxy clusters and groups has been on e of t he major achievements of the current 
generation of X-ray observatories (for a comprehensive review see |l45l | and references therein). Well-studied 
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. ^ L..g,.,^,[5MM, Hydra A [e.g., [HaiiTajTgjM^^ 

g.,[M[5l|92jM MI2Q0I, A 2052 [e.g. [m [2lHp , RBS 797 [e.g., IHlei, [lol llSi, A 133 [e.g 



examples of cavity clu sters are: Pers eus [e.g 



A 262 [e.g.,[2l|,[5i|, MS 0735+7421 [e.g., [To3,lT43^^ In-depth analyses of individual objects, which 



are now numerous in the literature, combined with studies of cavity samples [H, [Tt,, j64, 70, 72, 181j have 
enabled us to identify the global properties which are common among the cavities, thus shedding light of the 
feedback mechanism. The emerging picture is that bipolar outflows emanating from the BCG core inflate 
large bubbles while driving weak shocks, heat the ICM and induce a circulation of gas and metals on scales 
of several 100s kpc. Weak shocks have been observed as ripple-like features in the ICM in the deepest X-ray 
images of Perseus and A 2052 HH^. 

However, the differences between groups and clusters imply that the existing studies on cavities in clusters 
tell us little about how feedback operates in groups. With respect to rich galaxy clusters, the observation of 
cavities in galaxy groups and ellipticals is complicated by the lower X-ray surface brightness, which limit their 
detection in shallow X-ray images. On the other hand, there are several examples of AGN-ICM interaction 
just a few tens Mpc away from us which allow us to probe regions closer to the central black hole. In 
particul ar, low ni ass systems with cavities which now have deep Chandra images are: M 84 [91, 113] , NGC 
4636 [i,[iil[i6i], NGC 5044 ^,[11163,193, HCG 62 [loi HS, [HI Ell , NGC 5846 133]. Performing in- 
depth individual studies and sample studies of the lower-energy outbursts in these smaller systems is of major 
interest because the relationship between AGNs and hot gas can significantly influence galaxy evolution in 
the group environment, which is the locus of the majority of galaxies in the Universe j75| . Due to the 
shallower group potential, the AGN outburst can have a large impact on the intra-group medium in terms of 
altering the thermal history and spatial distribution of the intra-group mediu ni, as the mechanical output by 
radio AGN is of the same order of magnitude as the binding energy of gro ups [10l| . Such investigat ions have 



been undertaken o nly recently fo r individual objects (e.g., NGC 5813: [l83j |. AWM 4: [l67l Il69j ) and for 



group samples |67l . llOOl Ill7l l208j . and are rapidly improving our understanding of these systems. However, 
this observational effort is still awaiting detailed theoretical work in order to corroborate the observational 
findings. Recent detailed simulations indicate that groups are not simply a scaled-down version of clusters, 
as there may be remarkable differences between how AGN feedback operates in galaxy group and in galaxy 
cluster environments [l^j. In particular, AGN heating in groups seems to act through persistent, gentle 
injection of mechanical energy. On the other hand, in clusters there must be also the action of rare, powerful 
outbursts [94| . although more extensive theoretical work is required to reach firm conclusions. 
Examples of well-studied cavity systems in clusters and groups are shown in Figure [1] 



A. Cavity Heating 

The heating is thought to occur through the dissipation of the cavity enthalpy and through shocks driven 
by the AGN outburst. The energy required to create a cavity with pressure p and volume V is the sum of 
the pV work done by the jet to displace the X-ray emitting gas while it infiates the radio bubble, and the 
internal energy of the lobes, i.e. the enthalpy given by: 



E,,^ = H^ Eint +pV^ —^pV (23) 

where 7 is the ratio of the specific heats of the cavity content. If the lobes are dominated by the magnetic 
field, by non-relativistic gas, or by relativistic plasma, H can vary from 2pV and ApV. In particular, 
typically it is assumed that the internal composition of the cavity is dominated by relativistic plasma, 
therefore 7 = 4/3 and H = ApV. The product of pressure and volume can be estimated directly by X-ray 
observations through measurements of the cavity size and of the temperature and density of the surrounding 
ICM. A potential issue is represented by the uncertainties in the measurement of the cavity volume. The 
cavity size is usually estimated through a visual inspection of the X-ray images. This method is therefore 
dependent on the quality of the X-ray data, and also subject to the arbitrariness of the observer. The 
cavity size and geometry measured by different observers may vary significantly dependin g on the approach 
adopted, leading to differences between estimates of up to a factor of two in pV [e.g.. l47l . Il09l . Il68| . 

Systematic observational studies of samples of X-ray cavities show that their enthalpies measured from 
Eq. [23] lie between ~ 10^^ erg (in ellipticals, groups and poor clusters) and >10^^ erg (in rich clusters). On 
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FIG. 1: Top Left: The green contours outlining the 330 MHz radio emission from Lane et aL [l3ll ] are overlaid onto 
the 0.5-7.5 keV Chandra image of the galaxy cluster Hydra A (z=0.0538). The extended radio lobes fill a large-scale 
syste m of X-ray cavities and are surrounded by a "cocoon" shock. See also i]V Fl and i)VHI Adapted from Nulsen 
et al. [164| ]. Top Right: Very Large Array (VLA) radio contours overlaid onto the 0.5-7.0 Chandra image of the 
galaxy cluster RBS T9T (z=0.35). The subarsec resolution radio image shows the details of the innermost 4.8 GHz 
radio jets (blue contours), which clearly point in a north-south direction. Remarkably, these inner jets are almost 
perpendicular to the axis of the 1.4 GHz emission observed at 1" resolution (green co ntou rs) , which is elongated in 
the northeast-southwest direction filling the X-ray cavities. Adapted from Gitti et al. [l05| . Bottom Left: 0.3-2 keV 
Chandra image of the galaxy group NGC 5813 (z=0.0066) with 1.4 GHz VLA (blue) and 235 MHz Giant Metrewave 
Radio Telescope (green) radio contours overlaid. The image shows two pairs of cavities, plus an outer cavity to the 
northeast, two sharp edg es to the northwest and southeast, and bright rims around the pair of inner cavities. Adapted 
from Randall et al. [183 ] . Bottom Right: 235 MHz GMRT contours overlaid on the smoothed 0.5-2.0 keV Chandra 
image of the compact group HCG 62 (z=0.0137). The radio source shows a disturbed morphology with inner lobes 
clearly filling the well defi ned X -ray cavities, but with outer lobes having no associated X-ray cavities (see also W D|) . 
Adapted from Gitti et al. [lo|]. 



the other hand, simulations indicate that pV va ries with time during the cavity evolution and may be an 
inaccurate measure of the total energy released |l40l Il4ll | . Cavity power estimates within a factor o f tw o 
of the simulated values seem possible provided the inclination angle of the jets is known accurately [l49j . 
Bearing this caveat in mind, when divided by the cavity age, icav, the observational measurements give an 
estimate of the so called "cavity power", P^a y . S ince shocks are ve ry d ifficult to detect and are currently 
known only in a few systems (e.g.. Hydra A [ill], MS 0735+7421 [m^, HCG 62 [lo|, NGC 5813 flSSj), 
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for consistency the usual approach in sample studies is that of considering only the cavity power. Pcav thus 
provides a lower limit (and best-available gauge) to the true total mechanical power of the AGN, i.e. the jet 

power: Pjot >Pcav = ^^cav/^cav 

As proposed by Birzan et al. |16| . the cavity age can be estimated in three ways: (1) by assuming that 
the cavity rises the hot gas atmosphere at the sound speed Cg = jkT / finip: in this case the cavity reaches 
the projected distance R from the cluster center in the sound crossing time tg = R/cg = R/ ^kT / ^m-p-, (2) 
by assuming that the cavity is buoyant and move outwards at the terminal velocity Vt = \/2gV / SC , where 
g = GM^n/R"^ is the gravitational acceleration at the cavity position i?, V is the volume of the cavity, S 



is the cross-section of the cavity, and C — 0.75 is the drag coefficient [5l|: in this case the cavity age is the 



buoyancy-time tbuoy ~ R/ y^^gVjSC ; (3) by considering the time required for gas to refill the displaced 
volume of the cavity as it rises: in this case the cavity age is estimated as tref ~ 9, where r is the radius 
of the cavity. Typically the age estimates agree to within a factor of 2, with the buoyancy times lying in 
between the sound crossing time and the refill times. Most sample studies adopt the buoyancy time, which 
for typical values gives cavity ages of the order of few lO'' yr [e.g.. Il8l| . 

B. The Relationship between Jet Power and Lcooi 

Once a cavity is detected, it is relatively simple to estimate its power from the measurements of i?cav 
and tcav by applying Eq. 1231 The cavity power, Pcav, which is a measure of the energy injected into the 
hot gas by the AGN outburst, can then be compared directly with the gas luminosity inside the cooling 
radius, Lcooi, which represents the luminosity that must be compensated for by heating to prevent cooling. 
The gas luminosity inside the cooling radius is estimated as the bolometric X-ray luminosity derived from a 
deprojection spectral analysis. In Figure [5] (left panel) is shown a quantitative compar ison between Pcav — 



4py/ibuoy and Lcooi calculated for the extended sample discussed in O'SuUivan et al. 168j . who combined 
new data of 9 galaxy groups with the cluster sample of Birzan et al. [l7| and with the elliptical sample of 
Cavagnolo et al. [iF 



This plot follows those presented in Figure 2 of Birzan et al. [16| and in Figure 6 of Rafferty et al. [l8l[. As 
already noted by these authors, it is evident that the cavity power scales in proportion to the cooling X-ray 
luminosity, although with a big scatter. In general, it appears that the high mass (corresponding to high 
Lcooi) systems need an average of ApV per cavity to counter cooling. On the other hand, if we recalculate Pcav 
as IpT^/tbuoy all the points in the plot will shift down by a factor 4, and only the lower mass systems will still 
lie around the line Pcav = ^cooi- These systems require IpV per cavity to offset cooling at the present time. 
A few low mass systems will even still be above the equality line, thus indicating that the total mechanical 
power of the AGN far exceeds the radiative losses and their atmospheres are being heated. Although this 
extended sample is not a complete sample and therefore is not representative of the whole population of 
clusters, groups and ellpticals, it is interesting to consider the mean values of heating and cooling to see 
how they compare. We estimated a mean cooling power of 4.09 x 10^* erg s~^, and a mean cavity power of 
6.18 X 10^* erg s~^. In order to quantify properly the contribution of AGN feedback, over the system lifetime, 
in the energetics of cooling flow, it is important to determine the "duty-cycle" of AGN. Many studies have 
attempted such c alcu lation by adopting different approaches, e.g., by considerin g th e luminosity function 
of radio galaxies |l6lj . the fraction of cluste rs th at contain bubbles and cavities |70l-l72|, th e frequency of 
bubble required to produce sufficient heating [171^, the prevalence of radio- loud AGN [l5,[7i|. In particular, 
by considering the cavities as tracers of the feedback mechanism, i.e. by assuming that the feedback is 
active and efficient as long as the cavities are visible, we can correct the mean cavity power by the fraction 
of cool core clusters with cavities, estimated by Dunn and Fabian [t^I of the order of at least 70%. The 
ratio of mean cavity power to cooling flow power is thus very close to unity. The mean values for the whole 
sample are only indicative and do not reflect the different behaviour of groups and ellipticals with respect 
to clusters. In particular, such mean powers miss the point that in order to counter cooling the low mass 
systems require outbursts with relatively less total energy, lower powers and repeating more rapidly than 
high mass systems. This is supported by recent numerical simulations of galaxy groups which show that, 
in contrast to galaxy clusters, the AGN self-regulated feedback must act through a quasi-continuous gentle 
injection with subrelativistic outflows, rather than through rare and powerful episodes An attempt to 
produce more meaningful averages could be that of dividing the sample plotted in Figure [5] (left panel) in 
two subsets. In fact, by doing this we find that the ratio of mean cavity power to cooling power for the 
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FIG. 2: Left: Cavity power of the central AGN, Pcav, versus the X-ray luminosity of the ICM inside the cooling region, 
icooi- The cavity power is calculated assuming ApV of energy per cavity and the buoyancy timescale. Different 
symbols denotes systems in different samples presented in the l iterature: green triangles - Birzan et al. red 
squares - Cavagnolo et al. [4^, blue circles - O'SuUivan et al. [l68| . la uncertainties on cavity power are indicated by 
error bars (see 168i] for details). The diagonal line denotes Pcav = icooi- Credit O'Sullivan (private communication). 
Right: Cavity power of the central AGN, Pcav, versus integrated 10 MHz-10 GHz radio power, Lradio, for the systems 
in the sample of Birzan et al. (grey triangles) and the groups in the sample of O'Sullivan et a l. (black circles). 
The solid fit line indicates the regression fit to the data points calculated by O'Sullivan et al. [168| ]. The dotted line 
indicates the relation found by Birzan et al. [fTj. Adapted from O'Sullivan et al. ,168i] . 

groups and ellipticals is 7.94 (samples of Cavagnolo et al. [i^ and O'Sullivan et al. [l68l |). compared to a 
ratio of 1.51 calculated for clusters only (sample of Birzan et al. [i3|)- If the duty-cycle of low mass systems 
is the same as (or not lower than) high mass systems [70, 74], the relative ratio of heating to cooling appears 
to be a factor >5 higher in low mass systems. In other words, groups and ellipticals seem to have five times 
as much power available to counter cooling than rich clusters. 

A study of a complete, unbiased sample including both cool core and non-cool core systems is necessary to 
derive definite constraints on the balance between heating and cooling. However, it seems plausible that the 
time-averaged AGN feedback balances radiation losses of the ICM. Therefore the general picture emerging 
from the observed trend between X-ray luminosity and bubble mechanical luminosity, together with the 
existence of short central cooling time, is that the AGN is fueled by a cooling fiow that is itself regulated by 
feedback from the AGN. The basic idea of this AGN-cooling fiow scenario is that a self-regulated equilibrium 
may be achieved, in which the radiative losses from the thermal ICM are balanced by mechanical heating 
from the central AGN over the system lifetime. Although this scenario is no longer in doubt, is is still not 
clear how heating can act preserving at the same time the observed temperature gradient and the cool core 
[e.g-M- 

C. The Relationship between Jet Power and Radio Power 

Studies of cavity samples allow to derive the relationship between the mechanical power and radio emission 
of AGN jets and l obes . Su ch a relationship is of great interest because it helps to understand the physics 
of AGN jets [e.g., Il68l |226| |. and because it provides an estimate of the energy available from AGN based 
directly on the radio data [e.g., [l^, thus avoiding the problem of cavity detectability in shallow X-ray 
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images. Birzan et al. [17j studied a sample dominated by galaxy clusters and derived the relation between 
cavity power and 327 MHz radio power, as well as between cavity power and the integrated 10 MHz - 10 
GHz radio luminosity, extending to lower frequencies their previous work (l6| . While many of the observed 
X-ray cavities are filled with 1.4 GHz radio-emitting plasma, some are undetected at this high frequency 
and have been referred to as "ghost cavities". These may result from the aging of the relativistic particle 
population or may have been infiated by events which produced only particles of low energy. Examinations 
of radio images at multiple (and low) frequencies is particularly important as the progressive loss of particle 
energy causes higher frequency emission to fade fastest and the spectral index to steepen, so that evidence 
of a former AGN activity may be reflected only at low frequencies. The lack of 1.4 GHz radio emission is 
observed more frequently in groups than in clusters, therefore low- frequency radio observations are crucial 
for galaxy groups. 

Giacintucci et al. |l00j | selected a compilation of 18 galaxy groups, based on the presence of signs of 
interaction between the hot gas and the central AGN, and observed both by the Giant Metrewave Radio 
Telescope (GMRT) at frequencies < 610 MHz and by Chandra and/or XMM-Newton. These authors found 
that nine of thesegroups have cavities clearl y cor related with radio structures. By adding such systems to 
the Birzan et al. [l3| sample, O'SuUivan et al. [l68j | examined the relations between jet mechanical power and 
radio power in a combined sample which includes the groups having the most reliable radio measurements 
currently available. In particular, the integrated 10 MHz - 10 GHz radio luminosity estimated from the source 
spectral index is considered by these authors as a superior cavity power indicator compared to estimates 
at a single frequency, since it accounts for variations in spectral index between sources. Figure [5] (right 
panel) shows the relationship between cavity power, Pcav, and the i ntegrated radio luminosity, Lradio for 
the combined sample. The best fitting power-law relationship is [l68j |: 

log Pcav = 0.71 (±0.11) log Lradio + 2.54 (±0.21) (24) 

where Pcav and Lradio are in units of 10"*^ erg s^^. Sec O'SuUivan et al. |l68l | for a detailed discussion of 
this relation. 



D. Radio Lobe Composition 



When the radio source is filling the cavities, it is possible to compare the radio pressure of the relativistic 
plasma internal to the lobes with the X-ray pressure of the surrounding thermal gas. The pressure of the 
hot gas is measured from the density and temperature derived from the X-ray data as p ~ 2nckT. The total 
pressure in a radio lobe is the sum of the magnetic pressure, pb, and the total particle pressure, Ppart, and 
can be written as 

1 Ep,,t 1 (1 + k)E, 

Pradio = ±Ppart = + g "77^ " + Q 7T? '^^^> 

8tt 3/1/ Svr 3 j V 

where B is the magnetic field, k is the ratio of the energy in protons to that in electrons (E^), V is the 
volume of the radio lobe and / is the volume filling factor of the relativistic plasma. Using the expression for 
Ee given in Pacholczyk |l7l| . Eq. [25] determines the lobe pressure in terms of the magnetic field strength and 
the factor k/ f, once the volume V of the radio lobe is known. This calculation is usually performed under 
the widely adopted minimum energy conditions, in which the relativistic plasma is in equipartition with the 
magnetic field (Soq)- Further assumptions usually made in literature are / = 1 and fc = OorA:=l. A 
volume filling factor of 1 indicates that the lobes are empty of thermal gas, which is a reasonable hypothesis 
when they are observed to be spatially coincident with X-ray cavities. The assumption k = 1 implies that 
half of the energy in particles is in the form of non-radiating particles, as in an electron-proton jet, whereas 
k — would indicate an electron-positron jet. We stress that the uncertainties in the calculation of B^q 
and Pradio comc from the values of k and /, which are still largely unknown. Conversely, it is possible to 
constrain the ratio k/f by assuming pressure balance (see below). 

For historical reasons the frequency band adopted to calculate the standard equipartition field is i^i = 10 
MHz - j/2 = 100 GHz, i.e. roughly the frequency range observable with radio telescopes. From a physical 
point of view, the adoption of this frequency band in the calculation of the minimum energy is equivalent 
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to the assumption that only electrons emitting between 10 MHz - 100 GHz, i.e. with energy between 
7min oc (i^i/_Boq)^^^ and 7rnax cx: {^2/ B^q)^^'^ are present in the radio source. This approach neglects the 
contribution of the electrons emitting below 10 MHz and, as a more serious bias, in radio sources with 
different Scq selects different energy bands of the electron population because the energy of the electrons 
which emit synchrotron radiation at a given frequency depends on the magnetic field intensity [s^ . A more 
consistent approach is to calculate the minimum energy conditions, in which Boq does not depend on the 
emitted frequency band but directly on the low energy cut-off of the electron spectrum (typically assumed 
to be 7niin = 100). These so-called "revised" equipartition conditions select also the contribution to the 
energetics due to the low-energy electrons [st} . 

It is typically found that in cavity syst ems the X-ray pressure is more than one order of magnitude higher 
than the radio pressure [e.g., [l^ [s^, IH, Il09l |. It is also found that with revise d equipart ition the cavities 
are closer to pressure balance than they are with standard equipartition [e.g., Il09l Il68l |. Vice versa, by 
assuming the lobes are in pressure equilibrium with the ambient gas it is possible to constrain the particle 
content within the radio lobes [Ulea'j'eg'l. in particular, one can determine the ratio fcbai// that is required to 
achieve pressure balance under revised equipartition conditions. Several studies of the energetics and particle 
content of the radio lobes in cooling cores have found high values of fcbai//, up to several thousands (with 
standard equipartition) for active bubbles [e.g.,[l3, [z2|, suggesting that a large fraction of energy must be in 
non-radiating particles if / is close to unity. On the other hand, the pressure imbalance found in the lobes 
of FR-I radio galaxies in a sample of galaxy groups appears to be linked to the radio-source morphology, i.e. 
'plumed' sources typically have larger pressure deficits than 'bridged' sources where the jets are embedded in 
the lobes [5^. The authors interpret this result as evidence that plumed sources have a higher entrainment 
rate due to the larger fraction of the jet surface which is in direct contact with the external medium, leading 
to an increase in k/ f. Although the classification into bridged and plumed morphologies may not directly 
apply to radio sources at the center of cool core systems, typically having amorphous structures, this picture 
is consistent with the results of Dunn et al. [73] who argue that the large pressure imbalance observed in 
radio bubbles as those of the Perseus cluster is more likely to be due to entrainment rat her than a relativistic 
proton population. Recent studies show lobes having no associated X-ray cavities [e.g.. Il09lll68j |. Assuming 
their detection is not limited by the sensitivity of the current Chandra images, this suggests the possibility 
of mixing between ambient gas and radio plasma in the lobes. Therefore the /cbai// > values measured in 
such lobes is likely the results of entrainment of thermal gas through the hot gas atmosphere rather than an 
evidence of heavy jets ejected from the AGN. 



E. Radio Mini-Halos 



In some cases, the powerful radio galaxies at the center of cool core clusters are surrounded by diffuse 
radio emission on scales ~ 200 — 500 kpc having steep radio spectra (a > l^S^, oc v~°'). These radio 
sources, generally referred to as "radio mini-halos" , are synchrotron emission from GeV electrons diffusing 
through magnetic fields. Although the central radio galaxy is the obvious candidate for the injection 
of the population of relativistic electrons, mini-halos do appear quite different from the extended lobes 
maintained by AGN, therefore their radio emission proves that magnetic fields permeate the ICM and at 
the same time may be indicative of the presence of diffuse relativistic electrons. In particular, due to the 
fact that the radiative lifetime of radio-emitting electrons (~ 10* yr) is much shorter than any reasonable 
transport time over the cluster scale, the relativistic electrons responsible for the extended radio emission 
from mini-halos need be continuously re-energized by various mechanisms associated with turbulence in the 
ICM (reaccelerated primary electrons), or freshly injected on a cluster- wide scale (e.g. as a result of the 



decay of charged pions produced in hadronic collisions, secondary electrons). Gitti et al. 103j developed a 
theoretical model which accounts for the origin of radio mini-halos as related to electron re-acceleration by 
magnetohydrodynamic (MHD) turbulence, which is amplified by compression in the cool cores. In this model, 
the necessary energetics to power radio mini-halos is supplied by the cooling flow process itself, through the 
compressional work done on the ICM an d the frozen-in magnetic field. The successful application of this 
model to two cool core clusters (Perseus: jl03j and A 2626: |l04| ) has given support to a primary origin of 
the relativistic electrons radiating in radio mini-halos. 



Radio mini-halos are rare, with only about a dozen objects known so far. Gitti et al. |104| selected an 



initial sample of five mini-halo clusters based on the presence of both a cool core and a diffuse, amorphous 
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FIG. 3: Left: 1.4 GHz VLA radio contours (obtained by combining observations in A-, B-, and C- array configurations) 
overlaid onto the smoothed 0.5-2.0 keV Chandra X-ray image of the galaxy cluster RBS 797. The combined radio 
map has a resolution of 3", and is able to reveal the morphology of the central radio source, showing its elongation in 
the cavity direction (see top right panel of Figure[l]), without losing sensitivity at the larger scales. In particular, the 
extended radio emission is detected out to ~ 90 kpc. By subtracting the contribution of the central nuclear source, 
the residual flux density of the diffuse radio emission is ~ 11.5±0.6 mjy, indicating the likely presence of a radio 
mini-halo. Right: Integrated radio power at 1.4 GHz, [vP,^]-^ 4ghz' cooling flow power, Pc f = MkT / fim^, for the 
mini-halo clusters which have relevant X-ray and radio data available (data from [27l. IgsI. |99| . I104| . Il06|| i. 



radio emission with no direct association with the central radio source: Perseus 14l|, A 2626 [lOl flSSj. 
A 2142 [102], PKS 0745-191 [H, A 2390 In these clusters the size of the diffuse radio emission is 
comparable to that of the cooling region. These criteria are now typically adopted to identify mini-halos. 
However, the classification of a radio source as a mini-halo is not trivial: their detection is complicated by the 
fact that the diffuse, low surface brightness emission needs to be separated from the strong radio emission of 
the central radio galaxy. Furthermore, the criteria adopted to define mini-halos are somehow arbitrary (e.g., 
total size, morphology, presence of cool core) and some identifications are still controversial. This s aid, n ew 
detec tions of radio min i-halos ha ve re cently been claimed i n the galaxy clusters RX J1347.5— 1145 (l06j |. Z 
7160 nil, A 1835 [ml, A 2029 [m|, Ophiuchus [mllil^, RXC J1504.1-0248 pi], and RBS 797 (Fig. [3l 



left panel, see also [68 

Radio mini-halos are still poorly understood sources. Although secondary electron models have been 
proposed to explain the presence of their persistent, diffuse radio emission on large-scale in the ICM [e.g., 
Il26l Il74l | , a primary origin of radio mini-halos is now favored by recent statistical studies [i^l and by the 
observed trend between the radio power of mini-halos and the maximum power of cooling flows (see Figure 
131 right panel). This indicates a direct connection between cooling flows and radio mini-halos, i.e. the 
most powerful radio mini-halo s ar e associated with the most massive cooling flows, as expected in the 
framework of the Gitti et al. '103'] 's model. However, the origin of the turbulence necessary to trigger 
the electron reaccelcration is still debated. The signatures of minor dynamical activity have recently been 
detected in some mini-halo clusters, thus suggestir ig that a dditional or alternative turbulent energy for the 
reaccelcration may be provided by minor mergers pisl , Il06| and related gas sloshing mechanism in cool core 
clusters 143, 232]. Given the prevalence of mini-halos in clusters with X-ray cavities, another attractive 
possibility is that the turbulent energy is provided by a small fraction of the energy released by the bubbles 
rising from the central AGN (as suggested by Needless to say, a larger mini-halo sample as well as 

further theoretical investigations are necessary to reach a better understanding of this class of sources. 
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FIG. 4: Deep ~ 500 ks Chandra X-ray image (blue) and Very Large Array 330 MHz radio image (red) superposed 
with the Hubble Space Telescope visual image of the galaxy cluster MS 0735+7421. The giant X-ray cavities, filled 
with radio emission, are surrounded by a cocoon shock clearly visible in the Chandra image as an elliptical edge. The 
box is roughly 800 kpc by 800 kpc. 

F. Weak Shocks and Giant Cavities 

In addition to the cavity enthalpy, shocks driven by the AGN outburst may contain a large fraction of 
the energy released, thu s working to heat the ICM. Such shocks have been long predicted by numerical 
simulations [H, IH, Ill4| but are difhcult to detect since they are relatively weak (with Mach numbers 
~ 1 — 2) and are seen in projection against the cooler, brighter gas in cluster cores. We also note that to 
establish these surface brightness discontinuities as shocks one must measure an increase in temperature in 
the so-called "post-shock region" , as the ICM is heated by the passage of the shock . Usually the existing 
images are to o sh allow to rule out, e.g., the possibility that these features are cold front edges, due to gas 
slos hing [e.g., ll37l |. Besides a very few examples of strong shocks (e.g., Centaurus A, with Mach number ^8, 
[sgI . llSOj l). only recently elliptical surface brightness edges, consistent with arising from weak shock fronts 
driven by the cavities as they initially expanded, have become to emerge in deep Chandra expo sures of b right 
clusters and groups. Be autifu l examples of cocoon shocks are visible in the Hydra A cluster |l47l Il64| and 
in the NGC 5813 group ^8^, see left panels of Figure [T] 

Th e recent dis cover y of giant cavit ies and associated large-scale shocks in three clusters (MS 0735+7421 
[l47j |. Hercules A [l63j . Hydra A |164| ) has shown that AGN outbursts not only can affect the central regions, 
but also have an impact on cluster-wide scales possibly affecting the global properties of the ICM and the 
cluster scaling relations. In particular, the giant cavities discovered in the galaxy cluster MS 0735+7421 
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have a diameter of about 200 kpc each. The large volume of the cavities implies a huge cavity power: this 
large-scale outburst is the most powerful kno wn so far, releasing upward of 10^"'^ erg into the ICM and heating 
the gas beyond the cooling region [lOTl . Il47l | . The new, deep Chandra image has confirmed the presence of 
a weak (Mach number '^1.3) cocoon shock surrounding the cavity system (Figure H]). 

This new development may have significant consequences for several fundamental problems in astrophysics. 
As seen in SjTTl the observed relation between X-ray luminosity and gas temperature in clusters is steeper 
than expected if cluster growth were governed by gravity alone. This steepening is best explained by the 
addition of heat to the ICM and is therefore considered the main manifestation of nongravitational heating. 
The discovery of giant cavities has indicated that powerful AGN outbursts occurring at late times may 
contribute a significant fraction of the extra nongravitational energy. As mentioned above, this additional 
heating supplied by AGN could also induce the suppression of the gas cooling in massive galaxies required to 
explain the exponential turnover at the bright end of the luminosity function of galaxies [e.g.,[l3|- This would 
indicate a common solution for the two major heating problems associated with the ICM: those of cooling 
flow and galaxy formation. In the case of MS 0735-1-7421, th e dri ving energy of the shock as determined 
using a spherical hydrodynamic model is Es « 5.7 x 10^^ erg [147|. As estimated by Gitti et al. jl07l |. the 
AGN outburst in this cluster is heating the gas mass within 1 Mpc 7.7 x IO^^Mq) at the level of about 
1/4 keV per particle, and the heating level increases to ~ 0.6 keV per particle when considering the gas 

mass within 7'25oo: This is a substantial fraction of the 1-3 keV per particle of excess energy required to 

heat the cluster |228| . Only a few outbursts of this magnitude erupting over the life of a cluster would be 
required to heat it. By contrast, MS 0735-1-7421 is found to be a factor ~2 more luminous t han expected 
from its average temperature on the basi s of the observed L-T relation for galaxy clusters [l07l |. Based 
on the data presented in Gitti et al. [llOj |. we found a similar result for the giant cavity cluster Hydra A 
(Figure [5l left panel). Although caution should be taken in drawing general conclusions from the study 
of only a few objects, this indicates that flux limited samples of distant X-ray clusters may be biased in 
favor of detecting clusters with energetic AGN outbursts. We also note that powerful AGN outbursts may 
have a dramatic effect on the gas mass fraction measurements, due to an overestimate of the gas density 
[l07j |. The observed departure of MS 0735-1-7421 and Hydra A from the mean L-T relation is in apparent 
contradiction with the argument above that heat should steep the L — T relation, as also indicated by recent 
numerical simulations [e.g., 180]. However, we stress that the observed L — T relation is highly dependent 
on the definition of the characteristic temperature, i.e. for a fixed luminosity the position of each point 
in the plot may vary significantly depending on the choice of the method adopted to measure the average 
emission-weighted temperature for each cluster. Furthermore, the possibility of building a consistent L — T 
scaling relation from a sample of clusters relies on the capability to correct both the temperature and the 
luminosity measurements for the effects of the central cooling flow in a consistent manner for the whole 
sample. This may not be trivial as the physical conditions can vary significantly from case to case. For 
example, the commonly adopted method of excluding the central 70 kpc is found to have some drawback for 
giant cavity s yste ms as the cooling region and the effect of AGN feedback extend beyond this radius (see 
discussion in |l07l |). 

On the other hand, based on a study of XMM-Newton data, Gitti et al. (l07j | have shown that the energetic 
outburst in MS 0735-1-7421 does not cause dramatic instantaneous departures from the average properties 
of the ICM because it has not had a measurable impact on the large-scale te mper ature profile, which is 
in fact consistent within the scatter of the profiles observed in relaxed cluster j216l |. We recently found a 
similar result (Figure [SJ right panel) in the Hydra A cluster, although a sort of "plateau" standing below 
the typical profile in the range of radius ~ 0.05 — O.lrvir indicates the presence of cooler gas (see i jVHp . In 
general, these results suggest that there cannot have been many previous outbursts of high magnitude in 
these clusters, otherwise the total energy added to the ICM outside the cooling region should have had a 
marked effect. Studies of cavity samples found that the pre valen ce of outbursts as energetic as 10^^ erg is 3 
(namely, MS 0735-1-7421, Hercules A and Hydra A) over 30 [l8l| . If, as it appears from our in-depth studies 
of MS 0735-1-7421 and Hydra A, such powerful outbursts are rare in individual clusters, their occurrence in 
~ 10% of known cases hence requires that they occur a similar fraction of time in most cooling flow clusters. 
This picture is consistent with the model proposed by Nipoti and Binney |l6l| . in which the AGN activity 
is strongly variable with time and all systems occasionally experience powerful outbursts. 
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FIG. 5: Left: Bolometric X-ray luminosities corrected for the effect of cooling flow in the ce ntral ^ 70 kpc vs. 
emission- weighted temperatures derived excluding cooling flow components, from Markevitch |l36|. The dashed 
line is the best-fit power law to the sample: Lx = 6.35 • {kT/6keVf-^'^ x 10**ergs-\ The red triangle and the 
blue square represent the measurements from the XMM-Newton and Chandra data of the giant cavity cl uster s 
MS 0735-1-7421 and Hydra A, respectively, corrected consistently with the method adopted by Marke vitch [136| |. 
Right: Temperature profiles measured for a sample of relaxed clusters presented by Vikhlinin et al. (2161 ]. The 
temperatures are scaled to the cluster emission-weighted temperature excluding the central 70 kpc regions. The 
profiles for all clusters are projected and scaled in radial units of the virial radius r^ir, estimated from the relation 
rvir = 2.74Mpc^< Tx > /lOkeV [83|. Overlaid are the cooli ng fl ow corrected, scaled temperature profile s of the 
giant cavity clusters MS 0735-f7421 (red triangles, Gitti et al. [l07| and Hydra A (blue squares, Gitti et al. [llOl |). 



G. SMBH Growth 



A GN are powered by the release of gravitational binding energy from accretion onto massive black holes 
Il32l | . The matter that reaches the black hole converts it binding energy efficiently into AGN power as 
Pjet = eMc^, where e 0.1 — 0.4 depending on the spin of the black hole. Rapidly spinning black holes with 
spin parameters approaching unity are most efficient due to the smaller radius of the innermost stable circular 
orbit. The form of energy that is released depends on several factors including the accretion rate, the mass 
of the black hole, the structure of the accretion flow, and the spin of the black hole (see [l60j for a review). 
When black hole accretion approaches the Eddington rate, the binding energy is emitted thermally from 
an optically thick, geometrically thin accretion disk that is morphologically classified as a quasar or Seyfert 
galaxy. When the accretion rate drops below a few percent of the Eddington rate a radiatively inefficient 
AGN is formed (i.e., an ADAF), and the energy is released primarily in the form of mechanical energy 
associated with a radio jet. Accretion rates in bright AGN can be estimated using the radiation emitted 
from the nucleus that directly (e.g., UV or X-ray emission) or indirectly (eg., nebular or far-IR emission) 
trace the energetic output from the accretion disk [e.g.. Illi]. However, despite having mechanically powerful 



erg e 

radio AGN llal . brightest cluster galaxies (BCGs) rarely show strong X-ray and UV emission emerging from 



their nuclei 116] . implying that their accretion rates generally lie below a few percent of their Eddington 
rate. 



Rafferty et al. |181| estimated the accretion rates in a sample of BGGs in clusters with prominent X-ray 
cavities and found this to be the case. They estimated the accretion rates using the measured output power 
based on t he pV work done by the cavities divided by their buoyancy ages. Using this approach, Rafferty 
et al. |l8l| found that supermassive black holes centred in cool cores are growing at a rate of ^ 0.1 Mq yr~^. 
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In rare instances such as the powerful AGN in MS 0735+7421, [l07l Il47l Il48j |. the accretion rate exceeds 
1 Mq yr~^. Assuming black hole masses that are consistent with the values expected from their stellar 
luminosities and velocity dispersions, the accretion rates are consistent with being currently at most a few 
percent of Eddington. If AGN feedback in BCGs quenches cooling flows over the lifetimes of clusters, their 
black holes may be more massive than predicted by the A/bh ^ cr relation. 

H. Further Evidence for Mechanical Feedback 

As we have seen in WBl it is now widely accepted that the AGN in the cool cores can reheat the ICM. 
Although this is certainly its main impact, AGN feedback is likely to have other important effects on the 
ICM. We have recently investigated this point by performing an in-depth study of the galaxy cluster Hydra A, 
which harbors a well-known, large-s cale syste m of X-ray cavities embedded in a "cocoon" shock surrounding 
the central, powerful radio source ([146', '164"], see Figure [1] top left panel). By means of a detailed spectral 
analysis of the deep {'^ 240 ks) Chandra observations, we found indication of the presence of multiphase gas 
along soft filaments seen in the hardness ratio map (Figure [51 left panel) . Interestingly, such filaments are 
almost spatially coincident with the radio lobes of the powerful central radio source. The cooler gas has a 
significant impact on the radial temperature profile of the cluster, creating a sort of "plateau" which departs 
from the typical profile (blue squares in Figure[51 right panel). In fact, the scaled temperature profile of Hydra 
A measured after masking the filaments is found to agree well with the general shape of the temperature 
profi les observed for relaxed clusters, thus providing a confirmation that these filaments contain cool gas 
[llOj. By performing a spectral deprojection analysis of an absorbed 2-temperature component model, we 
found evidence that ^lO^^M© of low-entropy material has moved upward from the central 30 kpc to the 
observed current position of 75 — 150 kpc, likely due to some form of entrainment or dredge up by the rising 
lobes. Assuming that the mass of cool gas, which is ~60% of the total mass of gas remaining within 30 
kpc [ssj . was lifted out of the central cluster region by a continuous outflow or a series of burst s from the 
central AGN over the past 200 — 500 Myr (as it appears from the study of the cavity system, [227j ). it would 
amount to outflows of a few hundred Mq yr~^. There would thus be a development of gas circulation that 
can significantly reduce the net inflow of cooling gas, as initially discussed by David et al. [1^ and Nulsen 
et al. (l62j |. Therefore our results show that the AGN feedback in Hydra A is acting not only by directly 
re-heating the gas, but also by removing a substantial amount of potential fuel for the supermassive black 
hole (SMBH). This provides indications of mechanical AGN feedback acting through coUimated, massive 
outflows generated by jets or cavity dragging [e.g.. il78i |. 

The energy required to lift the cool gas gives a lower limit on the amount of AGN outburst energy deposited 
in the ICM. This value can be estimated by calculating the variation in gravitational potential energy during 
the lifting process. If we assume that the undisturbed ICM is approximately isothermal with sound speed 
Cg ~ 1000 km s ~^ and is in a hydrostatic configuration with density profile p(r), we can calculate this 
quantity as [l87l | 

^E^M^iJe£\ (26) 
7 \PfJ 

where Afcooi is the lifted mass, pi and p/ are the initial and final densities of the surrounding ICM, and 
7=5/3 is the ratio of specific heat capacities. From the density profile presented by David et al. (2001) we 
estimated that the energy required to lift the cool gas is >2.2 x 10®*^ erg. Th is value is comparable to the work 
required to inflate all of the cavities against the surrounding pressure |227l | and is ~25% of the total energy 
of the large-scale shock fl64j . We also note that a good fraction of the energy required to lift the low entropy 
gas will be thermalized when the gas falls back inward. Given the large energy required, uplift provides 
a signiflcant channel for the dissipation of outburst energy in Hydra A. There is a remarkable correlation 
between these low-entropy filaments and the meta l-rich filaments in the iron-abundance maps measured by 
Simionescu et al. [202] and Kirkpatrick et al. |l28l |. shown in Figure [6] (right panel). The emerging picture 
is that Hydra A's powerful radio source is able to lift cool, metal-rich gas from the central region and 
distribute it throughout the X-ray atmosphere of the cluster. A similar effect of extended m etal outfl ows 
in th e dir ection of the radio lobes and X-ray cavities is observed in other systems (e.g., M87 |200l . l223l [. A 
262: (l29j |. RBS 797 [1^). This is consistent with the results of the most current theoretical modeling of 
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FIG. 6: Left: High-contrast hardness ratio map of the galaxy cluster Hydra A obtained by dividing a 1.5-7.5 keV 
image by a 0.3-1.5 keV image. Regions in black are indicative of low temperature gas, indicating the presence of 
low-entropy filaments. Overlaid in green are the sectors used to study the spectral properties of the cool gas (located 
between radii ~ 70-150 kpc. We extracted the spectra in these sectors and compared two different spectral models: 
a single-temperature plasma "IT" model, and a "2T" model which includes a second thermal emission component. 
The F statistics for the spectral fitting improvement over the single-phase model indicate that the addition of a 
second thermal component is most significant in sectors SSE, NNW, NNE, and ENE, thus providing evid ence for 
the presence of multiphase gas in agreement with the hardness ratio map. Adapted from Gitti et al. [ll(j |. Right: 
Metallicity map showing the central 5' x 5' of Hydra A. Brighter regions represent a higher metallicity. The color scale 
of the metallicity Z (in solar units) is as follows: white: Z > 0.75, yellow: Z = 0.65 — 0.75, orange: Z = 0.55 — 0.65, 
red: Z = 0.45 - 0.55, blue: Z = 0.3 - 0.45, black: Z < 0.3. The 1.4 GHz radio emission is shown by the black 
cont ours. The green elliptical regions indicate the position of the inner cavities. Adapted from Kirkpatric k et al. 
[l28| . In both panels are overlaid the white contours outlining the 330 MHz radio emission from Lane et al. [l3l| . 



AGN feedback in massive cosmological systems, which predict the massive subrelativistic bipolar outflows 
and buoyant bubbles to produce a metal uplift along the jet axis (see W 



I. Numerical Simulations 



In the last decade the phenomenon of AGN feedback and its imp act on the ICM has been the subject of 
many theoretical investigations (see [186|, [l9l|, HI, [U, [HI, [Ha], [20|, [HI, [11, [11, [203, 34], 156], 



[95j , |9j] , to name a few) . These works focused either on several aspects of feedback physics and microphysics 
or with the global, long term evolution of the ICM. However, uncertainties are still large enough that the 
observations must guide researchers to select the relevant mechanisms at work in the feedback process. 

The body of different observational investigations set a number of key constraints on process of AGN 
feedback. The results presented in this paper strongly suggest that AGN feedback manifests itself as massive 
subrelativistic bipolar outflows which heat the ICM through weak shocks and form X-ray cavities, lift large 
masses of hot gas from the central region > 100 kpc and generate abundance asymmetries along the outflow 
direction. Processes such as AGN Compton heating or thermal conduction are unable to expl ain the collection 
of observations described above. Thus, although possibly relevant in some respect [e.g.. ll92| . they likely play 
a secondary role in local clusters and groups. 

Recent 3D hydr o simulations of outflow feedback ([93, ^B,', see, among others, also [H, 13, Il65l Il78l 

Il99l I207I I215I , l23l| ) have quantitatively verified that coUimated out flow s redu ce the gas cooling rate below 
the limits set by recent Chandra and XMM-Newton observations ( |l45l Il72j ] and references therein) for a 
timescale comparable with the cluster age. At the same time, this feedback mechanism generates ICM 
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density and temperature profiles which reasonably agree with those of typical cool core clusters. This is not 
a trivial result [SO, HH . The observable effects of the feedback from bipolar AGN outflows range from creation 
of X-ray cavities and large scale shocks to inducing entropy and metal abundance anisotropics, due to the 
lifting of central gas, relatively cool and metal rich, along the direction of the jet (see also Bruggen 2002; 
Roediger et al. 2007). These models, far from being complete and exhaustive (the origin of the outflows is 
essentially ad hoc), represent an interesting starting point for a more thorough understanding of the AGN 
feedback process. 

VI. CONCLUDING REMARKS 

In order to fully understand the growth and evolution of galaxies and their central black holes, the history 
of star formation, and the formation of large-scale structures, it is crucial to understand first the processes 
of cooling, heating and the dynamical evolution of the intra-cluster gas. In particular, the feedback from the 
central black holes has turned out to be an essential ingredient that must be taken into account in any model 
of galaxy evolution. The main manifestation of the action of AGN feedback is in galaxy clusters and groups. 
Their study, which is currently a very active line of research in extragalactic astrophysics, has allowed us to 
make significant progresses in this field. However, many details of the AGN feedback mechanism are still 
unclear. It is not well understood, for instance, how the heating distributes in space and time in order to 
drastically reduce gas cooling, preserving at the same time the central cool core. An even more puzzling 
issue is the process of black hole accretion and feedback energy generation. 

The last decade has represented a quantum leap in the quality of X-ray observations, thanks to the 
Chandra and XMM-Newton satellite telescopes. Recent results, discussed in this paper, have shaked up 
our understanding of the gas astrophysics in systems ranging from massive elliptical galaxies to rich galaxy 
clusters. They suggest that bipolar outfiows emerging from the BCG core inflate large bubbles, heat the 
ICM and induce a circulation of gas and metals on scales of several 100s kpc. 

The current generation of X-ray observatories is still working well and can be expected to continue doing 
so for few more years. Unfortunately, the prospects for the future of X-ray astronomy are not clear at the 
moment, and X-ray astronomers must rely on the good health on the existing, but aging, X-ray telescopes. 
As a result, the next few years represent a narrow time window to exploit the unique opportunity to observe 
deeply many additional objects, thus collecting crucial information on the cluster and group evolution. At the 
same time, current radio instrumentation is steadily improving both at the MHz and at the mm-wave ends 
of the spectrum, and is about to make a significant step forward with the next generation of observatories, 
such as Low Frequency Array (LOFAR) and Atacama Large Millimeter/submillimeter Array (ALMA). A 
common effort, from both the observational and theoretical side, will allow us to widen our knowledge on 
this fundamental problem which is central to the entire field of extragalactic astrophysics. 
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